Draft version December 14, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 10/09/06 



EVIDENCE OF STAR FORMATION IN LOCAL SO GALAXIES: SPITZER OBSERVATIONS OF THE SAURON 

SAMPLE 

Pasquale Temi 1-2 , Fabrizio Brighenti 3 ' 4 , William G. Mathews 3 

Draft version December 14, 2008 

ABSTRACT 

We discuss infrared Spitzer observations of early type galaxies in the SAURON sample at 24, 60 
and 170/im. When compared with 2MASS Ks band luminosities, lenticular (SO) galaxies exhibit a 
much wider range of mid to far-infrared luminosities then elliptical (E) galaxies. Mid and far-infrared 
emission from E galaxies is a combination of circumstellar or interstellar emission from local mass- 
losing red giant stars, dust buoyantly transported from the galactic cores into distant hot interstellar 
gas and dust accreted from the environment. The source of mid and far-IR emission in SO galaxies 
is quite different and is consistent with low levels of star formation, 0.02 - 0.2 M Q yr" 1 , in cold, 
dusty gaseous disks. The infrared 24/im-70/xm color is systematically lower for (mostly SO) galaxies 
with known molecular disks. Our observations support the conjecture that cold dusty gas in some SO 
galaxies is created by stellar mass loss at approximately the same rate that it is consumed by star 
formation, so the mass depletion of these disks by star formation will be slow. Unlike E galaxies, the 
infrared luminosities of SO galaxies correlate with both the mass of molecular gas and the stellar H/3 
spectral index, and all are related to the recent star formation rate. However, star formation rates 
estimated from the H/3 emission line luminosities Lap in SAURON SO galaxies are generally much 
smaller. Since Lap does not correlate with 24/im emission from dust heated by young stars, optical 
emission lines appear to be a poor indicator of star formation rates in SAURON SO galaxies. The 
absence of H/3 emission may be due to a relative absence of OB stars in the initial mass function or 
to dust absorption of H/3 emission lines. 

Subject headings: galaxies: elliptical and lenticular; galaxies: ISM; infrared: galaxies; infrared: ISM 



1. INTRODUCTION 

In the discussion below we extend our earlier analysis 
of Spitzer infrared observations of elliptical galaxies (E) 
(Temi, Brighenti Mathews, 2007a) to include lenticular 
galaxies (SO). The data set we consider here is further re- 
stricted to 31 early type E (19) and SO (12) galaxies that 
have been observed by both Spitzer and SAURON. The 
SAURON sample is generally representative and has re- 
ceived much attention from observers at all wavelengths, 
allowing detailed comparisons of Spitzer infrared obser- 
vations with many different galactic attributes: stellar 
populations, nature and extent of cold gas, dust content, 
etc. 

Our principal objective is to explore evidence of re- 
cent star formation in local SO galaxies at mid and far 
infrared wavelengths and to contrast this with the ap- 
parent lack of measurable current star formation in the 
sample E galaxies. Many SO galaxies in the Spitzer- 
SAURON sample have attributes consistent with recent 
star formation: significantly larger mid and far infrared 
luminosities, large masses of molecular gas, and younger 
stellar ages as inferred from optical absorption line spec- 
tra. Spitzer-SAURON SO galaxies having the largest far- 
infrared luminosities are known to contain rotationally 
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supported disks of molecular hydrogen with masses > 10 8 
Mq. Using the 24/i emission as a guide, stars are forming 
in many sample lenticular galaxies at very modest rates, 
0.02 - 0.2 Mq yr _1 . For some of these SO galaxies this 
level of star formation is comparable to the rate that gas 
is globally supplied from an old population of evolving 
stars, suggesting a quasi-steady state. In their optical 
spectra our sample SO galaxiess exhibit Balmer absorp- 
tion lines from young A-F stars, but the Balmer line 
emission from ionized gas associated with more massive 
O-B stars can be lower than expected. In addition, when 
star formation rates are estimated using standard empir- 
ical procedures, the rates determined from H/3 emission 
line luminosities in Spitzer-SAURON SO galaxies can be 
much less than those derived from 24/im luminosity. The 
cold, dusty molecular disks in some SO galaxies may orig- 
inate from stellar mass loss, may be a vestige of a more 
substantial disk during a previous incarnation as a more 
normal spiral galaxy, or may have been acquired by a 
merger event. 

2. THE SAURON-SPITZER SAMPLE 

In Table 1 we list 34 (non-dwarf) E and SO galax- 
ies in the SAURON sample that are also in the public 
Spitzer archive. This table contains additional details 
about each galaxy that will be discussed below. The 
SAURON sample is not complete but was chosen to in- 
clude nearby early type galaxies representing a limited 
range of absolute magnitude, ellipticity, environments, 
etc. as described in detail by de Zeeuw et al. (2002). The 
Spitzer-SAURON sample that we discuss here includes 
additional imponderable attributes imposed by infrared 
observers who contributed to the Spitzer archive. 



2 



Temi, Brighenti & Mathews 



3. OBSERVATIONS AND REDUCTION 

Far infrared data presented here were obtained with 
the Multiband Imager Photometer (MIPS) (Rieke et al. 
2004) on board the Spitzer Space Telescope (Werner et 
al. 2004) in three wavebands centered at 24, 70 and 
160 /zm. These data, collected from the Spitzer public 
archive, do not form an homogeneous and uniform data 
set in terms of image depth and observing mode. Table 1 
lists the program ID and PI name of the original Spitzer 
observing program from which sample galaxies have been 
selected. The reader is referred to these programs to ob- 
tain details on the observing modes and imaging strat- 
egy, as well as the on-source integration time for each 
target in the sample. As an example, Virgo Cluster 
galaxies recorded under the guaranteed time program 
(PID 69, PI G. Fazio) reach a relatively low sensitiv- 
ity of 0.5MJy/sr and l.lMJy/sr (la) at 70 and 160 fim, 
respectively, while other galaxies (i.e. PID 20171, PI P. 
Temi) have deeper maps at a sensitivity level of only 
0.12MJy/sr and 0.3MJy/sr for the same two wavebands. 
The SINGS data are recorded in MIPS scan mode, cov- 
ering a very large sky area (30' x 10' ), incorporating two 
separated passes at each source location. SINGS images 
correspond to maps with intermediate sensitivity. Apart 
from the SINGS observations, data have been acquired in 
MIPS photometry mode, allowing appropriate coverage 
of the sources and their extended emission. 

Reduction of Spitzer-S AURON data follows the same 
procedure used in Temi, Brighenti & Mathews (2007a). 
Here we briefly summarize the basic processes involved 
in the reduction. We started with the Basic Calibrated 
Data (BCD) products from the Spitzer Science pipeline 
(version 16.1) to construct mosaic images for all ob- 
jects. Final calibrated images have been produced us- 
ing the MOPEX (Mosaicking and Point-source Extrac- 
tion) package developed at the Spitzer Science Center 
(Makovoz et al. 2006). MOPEX includes all the func- 
tions and steps necessary to process BCD data into cor- 
rected images and co-add them into a mosaic. The ma- 
jor MOPEX pipeline used was the mosaic pipeline which 
consists of a number of individual modules to be run in 
sequence to properly perform the reduction. We refer the 
reader to the MOPEX webpage for a detailed description 
of each module. Since the data sets presented here have 
been acquired in different observing modes, the modules 
chosen in build-up the reduction flow and their parame- 
ter setup have been carefully selected to properly remove 
mode dependent artifacts in the final mosaic (outlier de- 
tection and median filtering). 

If foreground stars and background galaxies were 
present in the nal mosaiced images they were edited be- 
fore flux extraction was performed. These were identified 
by eye and cross-checked using surveys at other wave- 
lengths (Digital Sky Survey and 2MASS). Flux densi- 
ties were extracted from apertures that cover the entire 
optical disk (R25). Sky subtraction was performed by 
averaging values from multiple apertures placed around 
the target, avoiding any overlap with the faint extended 
emission from the galaxy. Statistical uncertainties re- 
lated to sky subtraction are usually less than 1% but can 
be appreciable (tens of percent) for faint sources. Ob- 
served infrared flux densities for each galaxy are listed in 
Table 2. Columns 2-4 contain the total flux F\A\ (mJy) 



at 24, 70 and 160 /im within each MIPS passband with 
respective widths AA = 5.3, 19 and 34 /im. 

Systematics in the MIPS calibration result in fluxes un- 
certain at the 10% level at 24 fim and 20% at 70 and 160 
/j,m. The uncertainties listed in Table 2 do not include 
the systematic uncertainties. Aperture corrections for 
extended sources were applied to the fluxes as described 
in the Spitzer Observer's manual. The corresponding 
MIPS luminosities are L\ = F\A\ ■ 4ttD 2 where dis- 
tances D are from Table 1. 

4. OBSERVATIONS 

4.1. Star Formation Rates from 24 Micron Emission 

Figure 1 shows the relationship between the 2MASS 
Ks band luminosities of galaxies in Table 1 and their 
MIPS passband luminosities at 24, 70 and 160 /im. The 
K s-band luminosities of E galaxies correlate tightly with 
L 2 4, but scatter dominates the correlations at 70 and 
160 fim. These results are similar to those discussed by 
Temi, Brighenti & Mathews (2007a) for a larger sam- 
ple of E galaxies. Since E galaxies typically have de 
Vaucouleurs surface brightness profiles at 24 /im (Temi, 
Brighenti & Mathews 2008), we conclude that most of the 
24 /zm emission originates in circumstellar dust around 
old red giant stars. However, Temi, Brighenti & Math- 
ews (2007a, b) argue that E galaxy emission at 70 and 
160/xm arises from colder, truly interstellar dust. But the 
stochastic variation that produces the scatter in the lower 
panels of Figure 1 is thought to be related to stochastic 
AGN energy releases from the central black hole that 
heats dusty gas in galactic cores which is then buoy- 
antly transported out 5-10 kpc into the hot gas atmo- 
spheres within these galaxies (Temi, Brighenti & Math- 
ews 2007a, b). Only for ellipticals with the lowest L 70 
and ii6o can the emission be explained by dusty mass 
loss from local stars. Rare E galaxies that are unusually 
luminous at 70 and 160/im probably have experienced re- 
cent mergers, but, as discussed below, there is no strong 
evidence for recent mergers in the SAURON sample E 
galaxies. Note that the vertical scatter of SO infrared 
luminosities in Figure 1 significantly exceeds that of the 
E galaxies at all wavelengths. This is one of the strik- 
ingly different properties of SO galaxies that we empha- 
size here. 

In Figure 2 the infrared luminosities are normalized 
by Lks (in solar unit) and the correlations at 24/im for 
E galaxies becomes much flatter. A least square fit to 
24/im emission from E galaxies in Figure 2 gives 

L 24 = 1 29-4±0-5 L 1.01 ± 0.04 crg (1) 

Evidently the nearly constant Log(L 2 4/ Lk s ) = 29.4 
for E galaxies represents the standard infrared SED of 
old stellar populations and their associated circumstel- 
lar dust. We show below that the departure of individ- 
ual E galaxies from this relation are uncorrelated with 
H/3 emission line luminosity. Therefore, it is unlikely 
that the small variations of cither infrared or optical line 
emission in SAURON E galaxies are related to apprecia- 
ble star formation or recent mergers with star-forming 
dwarf galaxies. 

While a few SO galaxies appear to have old stars with 
little or no current star formation, as in many E galaxies, 
others have significantly larger L^ILks which we inter- 
pret as evidence of a subpopulation of younger stars. We 
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imagine that interstellar dust in cold (possibly rotation- 
ally supported) gas is heated by massive, hot stars that 
recently formed from this gas. The large vertical scatter 
in the SO galaxies at all three wavelengths in Figures 1 
and 2 suggests that the interstellar dust in these galaxies 
has a much different origin and environment than far-IR 
emitting dust in E galaxies. 

In their study of star-forming galaxies Calzetti et al. 
(2007) propose a relation between the specific 24/im lu- 
minosity, (ALa)24 = A 2 4-La 24 , an< ^ tnc s ^ ar f° rma tion 
rate, 

SFR{M & yr- 1 ) = 1.24 x 10 - 38 [(AL A ) 24 (erg s -i)]0-8850^ 

(2) 

that can be used to find the SFR from global observations 
of (ALa)24- Values of (\L\)24 for our sample galaxies are 
plotted against Lk s in Figure 3. 

We shall assume that this SFR-(ALa)24 relation from 
Calzetti et al. holds for the SO galaxies in Figure 3, but a 
correction must be made for the relatively strong 24/im 
contribution from old stars in these early type lenticu- 
lar galaxies. Unlike the star-forming galaxies considered 
by Calzetti et al, where almost all of the 24/xm emis- 
sion evidently comes from interstellar dust irradiated by 
young hot stars, much or most of the 24/im emission 
from E and SO galaxies originates in circumstellar dust 
surrounding old, mass-losing red giant stars. To remove 
this background 24/im emission, we determine the SFR 
from the (ALa)24 emission in excess of that from the old 
stellar populations in our sample E galaxies collectively 
assumed to have SFR = 0. The correlation between the 
specific stellar 24/im luminosity for old E galaxies and 

Lks, 

Log(AiA)24,corr = (1 .01 ± 0.05)LogL Ks + 30. 1 ± 0.5, (3) 

is shown as a dashed line in Figure 3. Therefore to esti- 
mate the star formation rate in SO galaxies that lie sig- 
nificantly above this correlation in Figure 3, we enter the 
excess luminosity 

A(ALa) 2 4 = (ALa) 2 4 - (AL A )24, C orr (4) 

in the SFR-(ALa)24 relation above. The resulting SFRs 
for galaxies with positive ALog(ALA)24 are listed in the 
last column of Table 2. If SFR = is assumed for E 
galaxies, the mean scatter among the SFRs evaluate for 
elliptical galaxies in the upper section of Table 2, 0.014 
M ot iot yr -1 , suggests that the star formation rates for SO 
galaxies (in the bottom section) are uncertain to about 
±0.02 M odot yr- 1 . The SFRs for these local SO galaxies 
are all less than 1 M yr -1 . 

If the 70 and 160/im emission from SO galaxies is due 
to the same process that emits at 24/im, such as star 
formation, we would expect the 70 and 160/im luminosi- 
ties to correlate with 24/im residuals ALog(ALA) 2 4- This 
correlation is shown in Figure 4. Evidently, the same 
young stars in SO galaxies that heat nearby interstel- 
lar dust grains to emit at 24/im also heat more distant 
grains that emit at longer wavelengths. However, when 
the 24/im residuals for SO galaxies are plotted against 
70 and 160/im luminosities that are normalized with the 
Ks-band luminosities, as in Figure 5, the correlation is 
greatly improved. This figure suggests that Lfzr/Lks 
is likely to be an equally good measure of the star for- 
mation rate in SO galaxies as the 24/im residuals. Since 



Log(Ls/Lxs) (from Table 1) does not increase system- 
atically with 24/im residuals or SFR(24/im), we conclude 
that the low level of star formation in SO galaxies pro- 
duces little observable additional £>-band (or ifs-band) 
emission. Consequently, Lk s is a good measure of the 
mass of old stars and is insensitive to the small SFR in 
these SO galaxies. 

But Figure 5 informs us about a deeper relationship 
between star formation, nuclear molecular disks, and the 
mass of old stars in SO galaxies. Young, Bureau & Cap- 
pellari (2008) have studied compact CO-emitting molec- 
ular disks in four SAURON SO galaxies: NGC 3032, 
4150, 4459 and 4526. These disks typically have kpc- 
sized radii, total masses 0.5 — 5 x 10 8 Mq and mean sur- 
face densities 100 - 200 M Q pc~ 2 . In NGC 3032, the SO 
with the largest apparent star formation rate in Figure 
5, the molecular disk (disk radius: 1.5 kpc, mass: 5 x 10 8 
Mq) is counterrotating and must have been acquired by 
a galaxy merger in the past. The smaller molecular disk 
in NGC 4150 (disk radius: 0.5 kpc, mass: 0.55 x 10 8 M ) 
is kinematically oriented with a dusty disk that obscures 
optical light. However, NGC 4150 contains a counter- 
rotating core of younger stars having a spatial orienta- 
tion that differs curiously from that of the molecular gas 
and rotates in the opposite sense. This suggests that the 
molecular disk in NGC 4150 also resulted from one or 
more recent mergers. 

In the remaining two SO galaxies observed by Young, 
Bureau & Cappellari (2008) - NGC 4459 (disk radius: 
0.67 kpc, mass: 1.6 x 10 8 Mq) and NGC 4526 (disk 
radius: 1.1 kpc, mass: 5.7 x 10 8 M Q ) - the molecu- 
lar disks are well-aligned with the rotating system of 
old stars. These appear to be more normal SO galax- 
ies. NGC 4459 and 4526, together with NGC 3156 and 
3489, form a tight group of four SO galaxies in Figure 5 
having nearly the same Lfir/Lks and 24/im residuals, 
but with SFRs greater than normal elliptical galaxies. 
It is remarkable that among these four SO galaxies both 
Lfir and Lk s vary by about an order of magnitude 
yet they share the same location in Figure 5. Since the 
mass of old stars increases with Lk s and Lfir depends 
on the total mass of dust and star-forming cold gas, it 
is possible that the star-forming gas in these SO galax- 
ies is being resupplied by mass ejected from old evolv- 
ing stars. To explore this possibility, as Lk s varies over 
the range 10 10 — 10 11 (Figure 3), the rate that mass is 
ejected from old stars varies by about M* w 0.01 — 0.1 
Mq yr -1 where we assume (M*/Lk s ) = 0.60 (Bell et 
al. 2003) for an old stellar population with a Kroupa- 
like IMF and dM*/dt w 0.15M*/(10 n M Q ) M „ yr" 1 
(Mathews 1989). Within the uncertainties involved, we 
see that the mass supply from old red giant stars in these 
four star-forming SAURON SO galaxies (M* w 0.01 - 0.1 
Mq yr -1 ) is nearly equal to the SFR estimated from 
the 24/tm dust emission. Therefore it is possible that 
the internal production of gas in these four SO galaxies 
maintains rotationally supported disks of cold dusty gas 
that are forming young stars at the same rate, keeping 
the total mass of cold gas and dust rather constant with 
time. 

The two SO galaxies with the largest 24/tm residuals 
and L FIR /L Ks in Figure 5 - NGC 3032 and 4150 - have 
evidently acquired cold gas and dust in recent mergers, 
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and their current star formation rates exceed the rate 
that mass is supplied by stellar mass loss. As star for- 
mation proceeds in these galaxies, we speculate that they 
may evolve downward in Figure 5 (NGC 3032 NGC 
4150 — ► NGC 3156) until they reach the slowly-evolving, 
metastable situation described above for the four galax- 
ies including NGC 3156 where the SFR nearly balances 
the rate of stellar mass loss. This decreasing trend of 
star formation in cold nuclear disks may extend further 
to the only two elliptical galaxies in our sample known 
to contain molecular gas - NGC 2768 and 4278 - which 
also have the largest Lfir/ Lk.s among the E galaxies 
in Figure 5. When the stellar mass loss becomes com- 
parable to the SFR, progress along this downward trend 
in Figure 5 is expected to become slower. Finally, since 
the stellar mass loss rate in an old single stellar popu- 
lation decreases with time (oc t -1 - 3 ), we speculate that 
the oldest merger-free SO galaxies are those in Figure 5 
that have joined the cloud of elliptical galaxies in which 
the cold nuclear disks have no observable CO or far-IR 
emission. 

As seen in the infrared color-color plot in Figure 6 
the infrared color Log(L24/£7o) differs for those galax- 
ies in Figure 5 that are known to contain molecular gas. 
Since these mildly star-forming galaxies have the lowest 
Log(L 2 4/^7o) colors of the entire Spitzcr-SAURON sam- 
ple, this color could be used to identify other E and SO 
galaxies that contain similar molecular disks. 

4.2. Correlations of 24^m Emission with Hf3 Index and 
Molecular Mass 

From their analysis of SAURON galaxies McDermid et 
al. (2006) showed that the mean stellar age of SO galax- 
ies is systematically much less than that of E galaxies. 
We now show that this result is consistent with our ob- 
servations at MIPS wavelengths. Young stellar ages cor- 
respond to larger H/3 Lick (absorption) indices in optical 
spectra. In Figure 7 we show that the optical H/3 index in 
SO galaxies (evaluated within R e /8) correlates positively 
with bandpass luminosities at 24, 70 and 160 fim, sug- 
gesting that a substantial fraction of this infrared radi- 
ation comes from interstellar dust irradiated and heated 
by young stars. 

The correlations in Figure 7 show that the contribu- 
tion of young stars to the H/3 index in SO galaxies in- 
creases with Lfir/Lks since this latter ratio correlates 
with the 24/xm residuals (Fig. 5) and the SFR. We expect 
an even stronger correlation if H/3 is plotted against the 
24/^m residuals. This expectation is confirmed in the up- 
per panels of Figure 8 for H/3 indices determined within 
either R e /8 or R e of the galactic centers. SO galaxies 
with H/3 > 2 contain systematically younger stars. Fur- 
ther confirmation that ALog(ALA)24 is related to the star 
formation rate is provided by the lower panel of Figure 
8, showing that the mean stellar age from Table 1 de- 
creases with increasing residuals in Figure 3. The three 
E galaxies in Figure 8 with mean stellar ages less than 5 
Gyrs - NGC 2768 and 5831 - are somewhat anomalous 
in not having larger residuals ALog(AL,\)24- 

While observations of the HI gas mass in our Spitzer- 
SAURON sample are too sparse to consider here, 
more observations are available of their molecular mass 
(Combes et al. 2007). Figures 9 and 10 show that the 
molecular mass M# 2 in SO galaxies increases with all 



Spitzer infrared passband luminosities Lfir and also 
with Lfir/ 'Lks, although the scatter in Figure 10 is 
rather large (see also Combes et al. 2007, who found 
a similar correlation using IRAS data). A spurious unit 
slope correlation in Figure 9 could result in errors in the 
distances to the galaxies and this may account for the 
degrated correlation in Figure 10. Postivc correlations 
in these Figures would be expected if the emitting dust 
in SO galaxies is embedded inside molecular gas in which 
young stars are forming. The two elliptical galaxies con- 
taining molecular gas in these figures, NGC 2768 and 
4278, have small positive residuals in Figure 3, but are 
not otherwise exceptional. 

4.3. Disparity between SFR from 24^m and H(3 
Emission Luminosities 

One of the many advantages of SAURON integral field 
spectroscopy are improved measurements of the total lu- 
minosity of optical emission lines. In Figure 11 we plot 
the total luminosity in H/3 emission against far infrared 
luminosity in each MIPS passband normalized with the 
Ks band luminosity. E and SO galaxies appear to be- 
have in a different fashion in Figure 11. In our Spitzer- 
based sample of E galaxies we found positive correla- 
tions between the observed H/3 emission line flux and 
the bandpass flux Ffir at all three wavelengths, 24, 70 
and 160/im (Temi, Brighenti & Mathews 2007). How- 
ever, when the FIR luminosity is normalized with Lk s , 
as in Figure 11, the correlation between the H/3 emis- 
sion luminosity and L^ILks disappears for the E galax- 
ies. This indicates that the Ks — 24/^m color is part 
of the normal SED of an old stellar population and 
that the H/3 emission is uncorrelated with this color. 
At 70 and 160/rni, however, E galaxy Lup correlates 
with Lfir/ 'Lks- These correlations are consistent with 
the analysis of Temi, Brighenti & Mathews (2007a, b) in 
which both H/3 and FIR emission arise from buoyantly 
transported cooling regions widely dispersed in the hot 
interstellar gas of E galaxies. 

By contrast, the H/3 emission from SO galaxies in 
Figure 11 appears not to correlate with L^ /Lk s or 
Liqq/Lks and exhibits larger overall Lfir/ Lks than 
the E galaxies. Evidently H/3 emission from SO galax- 
ies in the SAURON sample is generally unrelated to far 
infrared emission from dust. The absence of this corre- 
lation for SO galaxies may seem surprising since massive 
young stars are the likely source of both optical line emis- 
sion and dust heating with subsequent far-IR emission, 
as in the Milky Way disk. Perhaps some of the H/3 emis- 
sion is absorbed by dust. 

In the top panel of Figure 12 we show this result in 
a different way by plotting the 24/xm residuals in Fig- 
ure 3 against the H/3 emission line luminosity. The ab- 
sence of any correlation for E galaxies indicates that there 
is no coordinated enhancement of emission from heated 
dust or optical emission lines that would be expected 
from star formation or from a recent merger with a star- 
forming galaxy containing gas and dust. This result is 
consistent with the finding of Temi, Brighenti & Math- 
ews (2005) that E galaxies with mean stellar ages (from 
Balmer indices) less than about 5 Gyrs often show no 
measurable enhancement at 24/zm. For E galaxies the 
H/3 emission is clearly unrelated to star formation. A 
low level of H/3 emission is always expected in E galaxies 
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from stellar mass loss into the hot interstellar gas that is 
ionized by hot post-AGB stars until it thermally merges 
into the hot interstellar phase. Larger, more stochastic 
contributions to optical line emission in E galaxies can 
result from cooling dusty hot gas that has been buoy- 
antly driven out of E galaxy cores by low energy AGN 
activity (Temi, Brighcnti, & Mathews 2007a, b). (Merg- 
ers as a source of strong E galaxy dust emission at 70 and 
160/im are also possible, particularly when there is addi- 
tional morphological or dynamical evidence of a merger 
event.) 

As before, SO galaxies in the top panel of Figure 12 
show no obvious correlation between Lup and 24/im 
residuals from Figure 3. The correlations of 24/im residu- 
als A(ALa)24, and therefore also SFR(24/xm), with both 
the stellar H/3 index (Figure 8) and the mass of molec- 
ular gas (lower panel of Figure 12) confirm that 24/im 
emission is likely to be a good measure of star formation 
in our sample SO galaxies. But the lack of a correlation 
between Lap and 24/zm residuals in SAURON SO galax- 
ies suggests that the standard observational method of 
determining the SFR from optical emission lines (e.g. 
Kennicutt 1998; Calzetti 2007) is inaccurate or inappro- 
priate for these galaxies which have relatively low levels 
of star formation. 

According to Calzetti (2007), the current star forma- 
tion rate can be determined from H/3 emission line lumi- 
nosities using 

SFR(Hf3) = 1.51 x 10- 41 L HI3 (ergs s" 1 ) (5) 

where we assume Ln a = 2.86Lnp since accurate H a line 
luminosities cannot be measured with SAURON. (While 
this Case B Ha/H(3 ratio is not expected to hold in the 
presence of dust absorption, it is unlikely to be seriously 
in error for escaping line radiation.) The SFR(H/3) in 
equation 5 is derived using a Kroupa-like IMF extend- 
ing to 120 Mq. Figure 13 shows that SFR(H/3) is much 
less than SFR(24/mi) for many of our SAURON sam- 
ple galaxies. For galaxies with the largest SFR(24/im) - 
NGC 3032, 4526 and 4459 - SFR(H^) is an order of mag- 
nitude smaller than SFR(24/im) . These are also the three 
galaxies in Table 1 with the largest masses of molecular 
gas, 1.6 - 5 x 10 8 M Q . 

Calzetti et al. (2007) also suggest a means of deter- 
mining the star formation rate in the presence of dust 
extinction that employs both Balmer line and 24/im 
luminosities, SFR(24/im,H/3). However, we find that 
SFR(24/mi,H/3) is very nearly identical to SFR(24/mi), 
i.e. the optical line emission is too low to affect the SFR 
from 24/im alone. 

If star formation occurs in intermittent bursts sepa- 
rated by more than the main sequence lifetimes of OB 
stars, ~ 10 7 yrs, the SFR from HII emission line lumi- 
nosities might often appear to be less than that from 
24/^m radiation. But it is very improbable that star for- 
mation in these massive (~ 10 8 M Q ) molecular disks 
would occur in globally coherent bursts. 

Alternatively, the relative absence of H/3 emission 
could be understood if massive stars are not forming in 
these SO galaxies; Calzetti (2007) notes that SFR(if/3) 
drops by ~5 if the maximum stellar mass is decreased 
from 120 to 30 Mq. However, the stellar initial mass 
function (IMF) may be radically altered in molecular 
disks in SO galaxies. Typical surface densities for the 



disks observed by Young, Bureau & Cappcllari (2008), 
100 - 200 M Q pc~ 2 , correspond to about 0.02 - 0.04 gm 
cm~ 2 which is less than the threshold density <~ 1 gm 
cm~ 2 required by Krumholz and McKee (2008) for mas- 
sive star formation. 

The column depths perpendicular to these molecu- 
lar disks correspond to hydrogen columns of about 10 22 
cm~ 2 for which the optical depth to dust absorption 
of Rf3 is about ten. Therefore it is conceivable that 
SFR(H/3) is underestimated simply because ionizing mas- 
sive stars form close to the central disk plane where line 
emission cannot escape. However, the 24/im residuals 
and therefore SFR(24/im) correlate with the stellar H/3 
index so at least some of the young stars with Balmer 
(absorption) lines must be visible. For this obscuration 
model to work it would be necessary for A and F stars to 
be dynamically excited away from the disk plane so their 
radiation can escape. Further insights about the miss- 
ing optical line emission in SO galaxies might be gained 
by observing emission lines further in the infrared, like 
for example Br7 (at 2.16/zm), which are less absorbed by 
dust. 

5. DISCUSSION 

The star formation rates in Table 2 found from 24/im 
emission using the SFR-(ALa)24 relation derived by 
Calzetti et al. (2007) suggests that the SFR in SAURON 
elliptical galaxies is 0.02 M© yr _1 and may in fact be 
zero. However, many SO galaxies in the SAURON sam- 
ple have SFRs that exceed 0.02 M yr -1 and extend up 
to almost 0.2 M yr" 1 for NGC 3032. 

Our results complement those of Young, Bendo &: 
Lucero (2008) who studied MIPS infrared emission from 
a sample of E and SO galaxies known to be very rich in 
molecular gas (M(H 2 ) « 1.5-64 x 10 8 M Q ). They found, 
as we do, evidence of enhanced 24/xm emission from many 
SO galaxies and showed that some SOs are measurably 
extended at this wavelength having sizes comparable to 
CO molecular images. This suggests that dust emitting 
at 24/zm is indeed associated with star formation in cold 
molecular disks. Three of the galaxies in the H 2 -rich 
sample of Young, Bendo & Lucero (2008) are also in the 
SAURON sample - NGC 3032, 4459 and 4526 - having 
M(H 2 ) = 2.6, 1.6 and 4.0 xlO 8 M Q respectively (Table 
1). 

Combes, Young & Bureau (2007) have detected CO 
emission from about 28% of E or SO galaxies in the 
SAURON sample. Of the ten galaxies with CO detec- 
tions shown in our Table 1, the majority (8) are SO galax- 
ies. They also find that CO-rich E and SO galaxies have 
stronger H/3 spectral indices. This is consistent with our 
finding that both M(H 2 ) and the H/3 index correlate with 
non-stellar 24/i emission in SO galaxies. 

Finally, the referee of our paper suggested that we also 
consider an analysis of the mid and far-infrared data 
in our sample using the new two-dimensional kinemat- 
ical classification of early-type galaxies proposed by Em- 
sellem et al. (2007) which is based on SAURON data. 
This new classification scheme is based on a single pa- 
rameter, X R = (R\V\)/{R{V 2 + cr 2 ) 1/2 ), where R is the 
distance from the galactic center in the image plane, V is 
the stellar radial velocity, a is the stellar velocity disper- 
sion, and the two-dimensional averages are taken within 
the effective radius R e . Galaxies with A_r > 0.1 are des- 
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ignated fast rotators and those having A# < 0.1 are slow 
rotators. Our Spitzer-SAXJTLO~N sample contains 9 slow 
and 22 fast rotators. We made alternative versions of 
Figures 1-12 using this new classification scheme and 
determined that all correlations and conclusions based on 
the figures presented here remain unchanged, although 
the E-like slow rotators are fewer in number and more 
sparsely distribued in the alternative figures. Neverthe- 
less, some plots are marginally different. For example in 
the Afl-versions of Figures 1-3, the slow rotators exhibit 
less scatter in the vertical direction, with fast rotators 
lying both above and below, but this may be a result of 
the smaller number of slow rotators (9) compared to E 
galaxies (19) in the ususal classification. Nevertheless, 
we do not include the alternative figures here because 
they contain no important new information. Further- 
more, we think it may be premature to abandon the fa- 
miliar morphological early-type classification scheme in 
terms of E and SO galaxies since the new Ar classification 
parameter is currently unavailable for the vast majority 
of early- type galaxies. 

6. CONCLUSIONS 

Perhaps the most important conclusion from these 
Spitzer observations of infrared emission from early type 
SAURON galaxies is that there are significant differences 
between SO and E galaxies in the relative content of dust 
and cold gas and star formation rates. Only a small sub- 
set of SO galaxies appear to share with E galaxies the 
same correlations between the 24/im luminosity, the H/3 
index, and molecular mass. Results from Sloan surveys 
of "early type" galaxies that combine SO and E galaxies 
(e.g. Kaviraj et al. 2007) may therefore be misleading 
if applied only to E galaxies. It is likely that the source 
of cold gas in E and SO galaxies is systematically quite 
different. Rotationally supported cold gas in SO galaxies 
may be a relic of their previous incarnation as late type 
spirals. 

The threshold for measurable star formation rates in 
the SAURON sample is about 0.02 M Q yr" 1 . All SO 
galaxies (except NGC 474, 1023, 7457 and 4382) exceed 
this value and all E galaxies (except NGC 4486 and pos- 
sibly NGC 2768) have lower SFRs. We have shown that 
the 24/um luminosity is a good measure of the SFR and 
that it correlates well with other related parameters such 
as the H/3 index, the mean stellar age determined from 
SAURON observations and the mass of molecular gas. 



Positive residuals from the E galaxy correlation be- 
tween AL\ at 24/im and Lks represent the emission from 
interstellar dust alone and are a measure of the galactic 
star formation rate. It is remarkable that these residuals 
for SO galaxies correlate so tightly with normalized far-IR 
luminosities Lqq/Lks and Luq/Lks- This suggests that 
the amount of dust, cold gas and old stars is rather finely 
tuned. The estimated mass loss rates from old stars in 
SO galaxies are very nearly the same as the star forma- 
tion rates derived from interstellar 24/im emission. This 
supports the notion that the mass of cold, dusty gas in 
some SAURON SO galaxies is in an approximate steady 
state - gas is lost by star formation at nearly the same 
rate that new gas is ejected from the dominant old stellar 
population. 

However, the H/3 emission line luminosity from SO 
galaxies appears to be very low and unrelated to the star 
formation rate. (This is also true for the global B — Ks 
colors of SO galaxies in our sample.) Star formation rates 
determined from the mid-IR luminosity (AL^)24 are of- 
ten about an order of magnitude larger than those found 
from the H/3 emission line luminosity, using the stan- 
dard SFR recipes. Perhaps star formation in cold gaseous 
disks in early type galaxies does not extend to massive 
OB stars because the disk column densities fall below the 
Krumholtz-McKee threshold of ~ 1 gm cm~ 2 . An alter- 
native interpretation is that OB stars do indeed form but 
are deeply embedded in dust that absorbs H/3 line emis- 
sion. If this is the explanation, radiation must still escape 
from slightly older A-F type stars since the star forma- 
tion rate based on 24/mi emission correlates with the 
Lick H/3 index in the spectra of young ~A stars. These 
non-ionizing stars may be less concentrated toward the 
central plane of the dusty molecular disks, allowing their 
radiation to escape. 
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hydrogen masses Mm are taken from Morganti et al. (2006). ^ Molecular gas masses Mh 2 are taken from Combes et al (2007). 9 Hf3 indexes 
(estimated within i? e /S) from Kuntschner et al. (2006). Ages are taken from McDcrmid et al. (2006). * Hj3 emission fluxes are taken from Sarzi 
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a The mean scatter among the SFRs for elliptical galaxies in the upper section of this table, 0.014 M ^ ot yr 1 , suggests that the (somewhat larger) 
star formation rates for SO galaxies in the lower section of the table are uncertain to ~ ±0.02 Mq yr - 1 . 
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Fig. 1. — SAURON galaxy luminosities in 24, 60 and 170/im Spitzer bandpasscs plotted against 2MASS Ks band luminosities. Detected 
E galaxies are filled circles and SO galaxies are open circles. Triangles represent upper limits. 
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Fig. 2. — Plot similar to Figure 1 but with T-*firI plotted against L^ s . The dashed line in the upper panel is a least squares fit 
for just the E galaxies (filled circles). The Ks - 24/mi color is nearly constant for old stellar populations in E galaxies regardless of Lk s - 
In general SO galaxies (open circles) have a much wider range of Lpm than E galaxies. 
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Fig. 3. — Plot of {\L\)2A against Ks-band luminosity Lk s , showing a tight correlation among E galaxies and uncorrected 24/j excesses 
for many SO galaxies. 
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Fig. 4. — Residuals ALog(AL^)24 in Figure 3 plotted against passband luminosities at 70 and 160/im. In this and subsequent figures 
Lfir represents either Lro or Lieo- 
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Fig. 5. — Residuals ALog(ALx)24 in Figure 3 plotted against passband luminosities at 70 and 160/iin normalized with Lks- Note that 
the correlation is much tighter than in Figure 4. 
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Fig. 6. — Infrared color-color plot for 24, 70 and 160/im emission from Spitzcr-SAURON galaxies. Labeled galaxies with known molecular 
disks lie systematically near the bottom of the plot with small Log(L24/^7o)- 
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Fig. 7. — Variation of T^firI ^Ks with stellar H/3 index (within R e /8). The Pearson correlation coefficients for the SO galaxies at 24, 
70, and 160/im are rso = 0.77, 0.61 and 0.46 respectively (discounting upper limits). 
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Fig. 8. — Upper panels: Residuals ALog(AL.\)24 in Figure 3 plotted against stellar H/3 index measured within R e /8 or R e . Symbols as 
in Figure 1. The Pearson correlation coefficients for the SO galaxies are rso = 0.77 for R e /8 and rso = 0.75 for R e . Lower panel: Residuals 
ALog(AI/A)24 in Figure 3 plotted against luminosity-weighted stellar age from Table f . 
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Fig. 9. — Variation of far-IR bandpass luminosities with mass of molecular gas. A positive correlation is seen at all three wavelengths, 
24, 60 and 170/im. Symbols as in Figure 1. 
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FlG. 10. — Plots similar to those in Figure 9 but with molecular mass plotted against far-IR luminosities normalized with Lk s - Symbols 
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Fig. 11. — H/3 emission line luminosities plotted against far-IR band pass luminosities normalized with Ly s . E galaxies (filled circles) 
show correlations at 60 and 170^tm, but SO data (open circles) show larger scatter at all three infrared wavelengths. The Pearson correlation 
coefficients at 24, 70 and 160/^ for SO and E galaxies are respectively: rso = 0.27, 0.51 & 0.32 and r E = 0.04, 0.78 and 0.72. 
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Fig. 12. — Upper panel: Residuals in Figure 3 ALog(AL^)24 plotted against H/3 luminosities. E galaxies are filled circles and SO galaxies 
are open. Lower panel: Residuals in Figure 3 plotted against mass of molecular gas. E galaxies are filled circles and SO galaxies are open. 
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Fig. 13. — Estimated star formation rates from 24/im emission plotted against star formation rates estimated from Hf3 emission line 
luminosities. Most SO galaxies (open circles) have SFR(H/3) << SFR(24/^m). Data is available for only 2 E galaxies (filled circles). 



